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Abstract 

We discuss the functional renormalization group and pion condensation in the presence of a finite 
isospin chemical potential /!/. We calculate the phase diagram as function of temperature T and 
fii. While the exact effective average action is invariant under certain gauge transformations, the 
effective action in the local-potential approximation is not. As a consequence, the critical chemical 
potential /ij for Bose-Einstcin condensation at T = is no longer equal to the mass of the condensing 
mode. We discuss possible solutions to this problem. 
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1. Introduction 

One of the goals of heavy-ion collisions is to create energy densities and temperatures high 
enough to create a quark-gluon plasma and to study part of the phase diagram of quantum chro- 
modynamics (QCD). The attempts to understand the phase structure of QCD, together with the 
sign problem at nonzero baryon chemical potential ^ts, have triggered interest in several QCD-likc 
theories. These include QCD at nonzero isospin density QCD with adjoint quarks, and two- 
color QCD . A common feature of these theories is that they are free of the sign problem. On one 
hand this admits their straightforward simulation using lattice Monte-Carlo techniques and serves 
as a check of these methods against model-independent predictions of chiral perturbation theory. 
On the other hand, it allows a direct test of various model calculations at high temperature and/or 
density where chiral perturbation theory is not applicable. Hence the study of QCD-like theo- 
ries contributes to our understanding of the physics of strongly-coupled gauge theories at nonzero 
temperature and density. 

In the present paper, we consider the simpler problem of pion condensation in the presence of 
a finite chemical potential fii for isospin using the linear sigma model. Lattice simulations [!, 0, HI 
suggest that there is a deconfinement transition of pions at high temperature and low density, and 
Bose-Einstein condensation of charged pions at high isospin density and low temperature. In fact, 
the deconfinement transition and the transition to a charged pion condensate seem to coincide. 
The deconfinement transition is found by measuring the Polyakov loop and the measurements 
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show a sharp increase (indicating deconfinement) at approximately the same temperature as the 
onset of pion condensation. The ground state energy and equation of state were studied on the 
lattice in la Jjfl- Pion condensation at finite /ij has also been studied using chiral perturbation 
theory [l|, [M ladder QCD [l(3|, the chiral quark model 11, 12\. the linear sigma mod' 
jil , Dyson-Schwinger equations |2(|, the Nambu Jona-Lasinio (NJL ) m odel 




MMEZl, and Polyakov-loop NJL (PNJL) models [H,[29|. In Refs. 




the effects of charge neutrality were also investigated. In contrast to the lattice simulations, the 
PNJL model suggests that deconfinement and the onset of Bosc-Einstein condensation are two 
different transitions. Finally, the possibility of a BEC-BCS crossover from a Bose condensate to 
superconducting state has been investigated [3(| . Note that this is not a phase transition since there 
is no global order parameter that distinguishes the crossover. Rather, it is a qualitative change from 
a pion condensate of tightly bound quarks to weakly bound Cooper pairs as the isospin density 
increases. 

The model calculations are typically of mean-field type and going beyond mean field is the 
next step in the study of these models. In order to do so we apply the functional renormalization 
group (FRG) [3l|, HI, M, [H 0, [H, |H, iH. The paper is based on Ref. 0. The functional 
renormalization group is a nonperturbative method, with a wide range of applicability. It is an 
alternative, but equivalent formulations to Wilson's ideas from the 1970s. It has been used success- 
fully to calculate critical exponents for phase transitions 14 1| and to map out the phase diagram 
as function of temperature and baryon chemical potential |42. 43l 44, 45 as well as the calcula- 
tions of thermodynamic functions [461 147| and momentum-dependent correlation functions at high 
temperature jijj . 



2. Flow equation and derivative expansion 

One implementation of the renormalization group ideas is based on the effective average action 
Tk [<p\ as proposed by Wetterich (3l| . The effective average action is a functional of a set of fields 
denoted by 4> and satisfies a flow equation which is an intcgro-diffcrcntial equation. The subscript 
k indicates that all momentum modes q between the UV cutoff of the theory, A, and k have been 
integrated out. When k = A, no momentum modes have been integrated out and the effective 
action is equal to the classical action of the theory, i.e. Ta[</>] = S[<f>]. When k = 0, all modes 
have been integrated out and Tq [4>] is equal to the full quantum effective action. All quantum and 
thermal fluctuations from k = up to k = A have then been included. The exact flow equation for 
T k [4>] is 



1. 



d k T k {cp] = -Tv[d k R k (q) [r^ + Rk(q)] qt _ q \ . (!) 

where the superscript n on rj^[</>] means the n'th functional derivative of T k [(j)] and the trace is 
over the spacetime momenta q, and indices of the inverse propagator matrix. The function R k (q) 
is a regulator and is introduced in order to implement the renormalization group ideas: R k (q) is 
large for q < k and small for q > k whenever < k < A, and R\(q) = °o. These properties ensure 
that the modes below k are heavy and decouple, and only the modes between k and the UV cutoff 
A are light and integrated out. The diagrammatic representation of the flow equation is shown in 

Fig.rj 

The choice of regulator has been discussed extensively in the literature m, sum [m [ii . we 

employ the regulator which is given by 



2 



dk^k — ^ 



Figure 1: Diagrammatic representation of the exact flow equation for the effective action [<j>]. The line denotes 
the exact field-dependent propagator and the circle denotes an insertion of the regulator Rk(q)- 



Rk(q) 



(fc 2 - q 2 )6(fc 2 - q 2 ) 



(2) 



This regulator was introduced in Ref. |53[ (see also [5l|). and is particularly convenient in practical 
calculations since one can carry out the integral over three-momenta exactly and this turns the flow 
equation into a partial differential equation. 

The starting point of our discussion is the classical action S [(/>]. In the present case we consider 
N complex scalar fields with a quartic self- interaction and S[<p] reads 



S[<f>] 



where the classical potential is 




$| +$i 



V(4>) 



(3) 



V(4>) = 



2N \ i ii 



(4) 



and where i = 1,2, 3...-ZV. The complex fields $j are written in terms of real fields, 4>j = (4>2i-i + 
*</ > 2i)/v / 2- We identify tfii with a, (02±i<fe)/\/2 with tt^, and 04, 4>5, <p2N with 2^ — 3 neutral tt° . 
The parameter H breaks the 0(2N) symmetry explicitly to 0(2N — 1) and gives rise to nonzero 
pion masses. The parameter m 2 is negative in the remainder to ensure spontaneous symmetry 
breaking in the vacuum for H = 0. 

A chemical potential //* that corresponds to the conserved charge of the (7(l)-symmetry — > 
e lc *$i is introduced by the replacement do&i — > (do — The chemical potential acts as the 

zeroth component of an Abelian gauge field. The classical action is therefore invariant under the 
gauge transformation 



Hi + idoa . 



(5) 
(6) 



The introduction of the isospin chemical potential fii reduces the symmetry to 0(2) x 0(2N — 2) 
if H = 0, otherwise to 0(2) x 0(2N — 3). Chiral symmetry is broken in the vacuum, either 
spontaneously or explicitly by a vacuum expectation value (f>o for <fii. We also allow for a pion 
condensate by a nonzero value po of fa- A pion condensate breaks the 0(2)-symmetry. 
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After having introduced the isospin chemical potential pi, the classical potential is 



1 aW>o + Po) 



t -9a i 
2 PlPo 8N^° Po ' 



2\2 



Ufa , 



(7) 



where the subscript A is a reminder that the classical potential equals the boundary condition for 
the effective potential £4(0), i. e. Ua(<P) — V{4>) ( see below). 

The exact flow equation ([!} cannot not be solved so one must resort to approximations. The 
derivative expansion is a commonly used approximation and reads 



dr 



Z\ (V^) 1 ■ (V$i) + Z\ [(do - in) Po - pi) $i] 
+Zl (V$ 2 ) f • (V$ 2 ) + ... + E/feO) + 



(8) 



where Z\ are wavefunction renormalization constants and Uk{4>) is the scale-dependent effective 
potential that depends on the classical fields that are collectively denoted by (j>. Note that the 
Z^s are in principle different since Lorcntz invariance is broken due to both finite tempera- 
ture and finite density. The ellipsis indicate higher-order operators including derivatives such as 
${<t>i [(V$i)^ ■ (V<f>i)] . In the local-potential approximation, we set Z l k = 1 and omit such higher- 
order operators. In this approximation, the flow equation is an integro-diffcrential equation for 

The flow equation for the effective potential Uk{4>) is derived in the Appendix and reads 



d k U k {4>) 



Jd~l) 



IsfP,. 



+ (2N - 3) 



l + 2n(u 3 ) 

2^3 



(9) 



where d— 1 is the dimension of space, Vd = [2 d+1 7r <i / 2 r(i)] -:L , n(x) = l/{eJ 3x — I) is the Bosc-Einstein 
distribution function, Pi and Ri arc poles and residues defined in the Appendix. Moreover, Uk(4>) 
is a function of p\ — \4>\ and p 2 = \p\- We then introduce the shorthand notation diUk(4>) = 
dU k (<f>)/d Pl , d 2 U k (cf>) 



dU k (<f>)/dp 2 , and uj 3 = y/k 2 + diU k . 



3. Numerical results and discussion 

In this section, we solve the flow equation ([9J numerically using a third-order Runge-Kutta 
method and discuss the resulting phase diagram. The boundary potential U\(<j>) is given by Eq. ([7]), 
where we in remainder of this section set N = 2. The parameters m\ and g\ are determined so 
that we reproduce the pion mass m^, the sigma mass m a , and the pion decay constant f„ in the 
vacuum, i.e. for T = p,j = 0. In the vacuum we have 4>q — /„■ and po = 0. The A:-dcpcndcnt masses 



1 Note that in the vacuum and for H = 0, the potential depends only on combination + p\- One can therefore 
use the symmetry to rotate away po- This is in accord with the Vafa-Witten theorem [o-ll that states that parity 
cannot be spontaneously broken in the vacuum. 

2 Somctimes one further assumes that U^{<f>) takes a polynomial form. The partial differential equation (J9j is then 
turned into a set of coupled differential equations for the couplings. 



4 



m n h and m ak can be related to the effective potential at the fc-dependent minimum f n k as follows 



l TT,k 



diU k ((f>) 
diU k {<j>) 



(10) 



*>0 = /x,fc,PO=0 



ml k + <fi%U k (<f>) 



f>a=fT,,k,Po=0 



(11) 



where Uk(4>) = Uk(4>) + Hcf>o, i.e. the potential without the explicit symmetry breaking term. In 
particular, the physical pion and sigma masses, i.e. the masses for k = are given by 



t>0 = f7r,PO=0 



(12) 
(13) 



0O=/tt,PO=O 



We can now use Eq. (fTTj) together with <^o,fc = fn,k to determine the two parameters m\ and g\. 
The fact that = implies that this parameter is already known and the term —H(f>Q can be 

added to the potential after we have determined m\ and g\. 
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Potential at k = for fii = MeV. 

■ Potential at k = for fii = 30 MeV. 

■ Potential at k = A (/t/ = MeV). 
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Figure 2: The real part of Uk(<f>o,PO = 0)/A 4 as a function of 0q/2A 2 at T = in the chiral limit. The potential 
is plotted for k = A, and for k = in the case of /i/ = and = 30 MeV, where we have neglected the constant 
added to the potential by renormalization. We have used m<r = 400 MeV and f n = 93 MeV in order to tune the bare 
parameters of the potential at the cutoff A = \fZm a . The minima of the renormalizcd potentials have been marked 
with open circles as a guide to the eye. 
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In Fig. [2] we show the real partH of the scaled effective potential as function of </>q/2A 2 in the 
chiral limit and for T = at /xj = and /Ltj = 30 MeV. 

We use the experimental value f v = 93 MeV and m a = 400 MeV to tune the bare parameters 
in this case. We chose the cutoff A 2 = 5m 2 . This gives the values m\ = — 0.96A 2 and g\ = 37.56. 
We find that the inclusion of a finite chemical potential has the effect of increasing the po value 
for which Uk(4>) takes its minimum. We also see that including the quantum fluctuations from 
the bosons reduces the amount of symmetry breaking. For fermions, the effect is opposite due a 
sign change of the vacuum term in the flow equation. Moreover, if the symmetry is restored for 
some values of T and at k = 0, there exists a scale k = fcsSB such T still displays symmetry 
breaking when k > &ssb- 



> 

QJ 



60 



40 



20 









\ \ / <• 
\ r * 

Y v.,/ • 




'•X N --. 

/ V . 




i \ 

-' " 


1 1 


i ■ 

i 

/ i 



50 



100 150 

VI (MeV) 



200 



Figure 3: The chiral condensate and the charged pion condensate as a function of isospin chemical potential at the 
physical point at T = 0. Solid line is the local-potential approximation with N = 2, dashed line is the large- N limit 
and the dotted line is the mean-field approximation. We have used the values m CT = 600 MeV, = 140 MeV, and 
= 93 MeV, and set A = 600 MeV. 

In Fig. [3j we show the chiral condensate and pion condensate as function of [ii at zero temper- 
ature. We have used the parameters = 93 MeV, m a = 600 MeV, = 140 MeV, and A = 600 
MeV. The solid line is the local-potential approximation with N = 2, the dashed line is the large- TV 
result, and the dotted line is the mean- field result. In the chiral limit the charged pion condensate 
is the only condensate present at nonzero chemical potential when T = 0. If we go to the physical 
point however, we see that a chiral condensate appears in the direction of 4>q. For small chemical 
potentials, we find that the chiral condensate is the only existing condensate, that is, the minimum 
of Uk(4>) lies on the (/>o-axis and no charged pion condensation takes place. As we increase the 
chemical potential however, the minimum eventually moves away from the c/)o-axis and a charged 
pion condensate is formed, i.e. the chiral condensate is (partly) rotated into a condensate of charged 
pions. At T = we find that this happens when the chemical potential is about the vacuum pion 



3 Note that the effective potential in the broken phase has an imaginary part to the left of the minimum. 
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mass. This is in agreement with the results of e.g. [l(J, 12, 16, 17, 26| (but see the discussion below) 



K* 150 



Physical point 

Chiral limit 
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Figure 4: Phase diagram for the charged pion condensate in the chiral limit and for the physical point. We have 
chosen physical values = 93 MeV, = 600 MeV and m n = 140 MeV (m-^ = in the chiral limit) in order to 
tune the bare parameters at T = 0. We have also set A = 600 MeV. The symmetric phase is above the lines, while 
the broken phase is below the lines. 



In Fig. |4l we plot the phase diagram as function of temperature and isospin chemical potential 
in the chiral limit (dashed line) and at the physical point (solid line) . We use the same parameters 



me paramct 

as before. The phase diagram is in qualitative agreement with previous results [12|, 3,113, 0|- In 
particular, the phase transition to a condensate of pions is of second order and this is in accord with 
universality-class arguments. Note, however, that the critical chemical potential for T = and at 
the physical point is different from the pion mass m T and this is in disagreement with exact results. 
In fact, it can be shown that the critical chemical potential for Bose condensation at T = is equal 
to the mass of the condensing mode. Using the bare parameters we found at T = /i/ = 0, we 
see from Fig. 2] that we obtain a slightly different value for the critical chemical potential, namely 
/ij = 135 MeV. The reason we do not get the pion mass is that we are working in the local-potential 
approximation, where wave-function renormalization is neglected. Recall that to include an isospin 
chemical potential we should let <9o<& — > (do — /Uj)$, where we view the Lagrangian as a function 
of complex field variables. This term remains unrenormalized in the local-potential approximation 
and so the effective action is not gauge invariant under the transformations j5]) - j6]). This is in 
contrast to the full quantum effective action at T = [HBT ]. 

A simple solution to this problem is to first solve the flow equation for the effective potential 
Uk{4>) in the absence of the isospin chemical potential and then gauge the effective action in the 
local-potential approximation by the substitution do& (do ~ The effective potential at 
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scale k may then be written as 



U k {4>) = U k {4>)-\lJL 2 I pl-H(t> Q , (14) 

where U k {4>) denotes the 0(4)-symmetric potential, i.e. the potential satisfying Eq. ((9j) with fii = 0, 
that is Eq. (f3~4"|) in the Appendix. We call this the 0(4)-approximation. In order to get a charged 
pion condensate at k = 0, one needs a negative partial derivative in the po direction at the minimum 
of U k (4>). When [ij is just large enough for this to happen, we obtain 

= d 2 U k=0 (4>) 

= d 2 U k (cl>)-ri. (15) 

Using the 0(4)-invariance of U k (4>) and combining Eq. (IT5l) with Eq. (fT2|) . this yields a critical 
chemical potential 

Hj = m v . (16) 

While the effective potential U k (4>) defined by Eq. (|14[) is in agreement with exact results, its 
dependence on the isospin chemical potential is of course trivial. 

The large- N limit provides us with another flow equation for U k {<p) which is independent of the 
chemical potential. It is given by Eq. (|35[) in the Appendix: 

d k U k = ^ Vd - 1 [1 + 2»(W3)] . (17) 

(d - l)L0 3 

Thus using the same arguments as above, we obtain the critical chemical potential for pion con- 
densation which is equal to the pion mass. 

In Fig. [5l we show the phase diagram obtained by solving the flow equation for U k {4>) in the 
0(4)-approximation and in the large-iV limit. For comparison, we also show the phase diagram 
obtained by solving Eq. © and was shown in Fig. 2] albeit with a different cutoff A. In this plot, 
wc used A = 1500 MeV (Note that the dotted lines arc from Fig. H with A = 600 MeV). We would 
like to point out that both the 0(4) approximation and the large- ./V approximation yield a second- 
order transition and adding the term — \[iip\ does not change this. Hence all three approximations 
predict the correct order of the transition. 

In closing, we remark that it would be interesting to include the effects of wavefunction renor- 
malization. This would lead to a set of coupled equations for Z\ and U k ((j)). Hopefully it would close 
the gap between the value /ij = 135 MeV obtained in this work and the exact result pj = m OT = 140 
MeV and so perhaps provide a practical solution to the problems of the local-potential approxi- 
mation. Similarly, it would be interesting to apply the method to the chiral quark model. Since 
this model has quark degrees of freedom, one can simultaneously study the effects of a baryon 
and an isospin chemical potential. In particular, one can investigate the BEC-BCS transition as 
one increases pi. For small values of /x/, but larger than m w , there is a condensate of pions. For 
asymptotically large values of /Uj, and therefore large densities, the system is weakly coupled due to 
asymptotic freedom. In that case, the weak attraction gives rise to Cooper pairs and a BCS state. 
The order parameter has the same quantum numbers as the BEC state and so the transition is a 
crossover [l[. 

The BEC-BCS transition has been discussed at length in the functional renormalization group 



framework for ultracold atomic gases [56l . 1571 . l58| . The crossover in these nonrelativistic systems is 
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Figure 5: The phase diagram for the 0(4)-approximation (solid line) and the large- N limit (dashed line), where we 
have set m v = 140 MeV, f n = 93 MeV, m a = 600 MeV, and the cutoff A = 1500 MeV to determine the parameters 
and g\. We have also plotted the previous obtained phase diagram for comparison (dotted line). 



realized by changing the s-wave scattering length. This is done by varying an external magnetic 
field B and exploiting the existence of a Feshbach resonance at Bq = 834 Gauss. In this way one 
can basically change the coupling constant from minus infinity to plus infinity. A BEC state is 
obtained for positive scattering lengths and a BCS state for negative scattering lengths, i. c. on 
each side of the resonance. 
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4. Appendix 

In this Appendix, we derive the flow equation ([§]). The starting point is the exact flow equa- 
tion mj 



-Tr 



d k R k , q [T^ + R k {q)] \ 



(18) 



The effective potential Uk(<f>) is defined by evaluating the effective action r &[</>] for space-time 
independent values of the fields, i. e. 



Uk(4>uni) 



1 

VT 



(19) 
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where VT is the spacetime volume of the system and </> un i is constant. Taking the derivative of 
Eq. (Ti"9| with respect to k, we obtain the RG-equation for the effective potential 



dkU k ((j) uni ) 



1 



2VT 



:Tr 



dkRk, q [r 



(20) 



In the remainder we drop the subscript uni. With T k [(f)] given as in Eq. (|18[) . we find the Fourier 

.(2 
k 

" d 2 U k 



( 2) 

transform of TV [<f>] to be 



,(2) 

k,i,j,q,q' 



(2ir) d 6(q + q') 



+ Sijq - 2fj,rq (Sj2Si3 - 5 t2 S j3 ) 



(21) 



where d denotes the spacetime dimension and i,j are the indices of the propagator. Note that by 

(2) 

the presence of the chemical potential, T k t . , is not a diagonal matrix. Also note that we cannot 
rotate the field vector </> = (</>i, .., 4>2n) to point in one specific direction as we do not have 0(2N) 
symmetry. This has been explicitly broken to a 0(2) x 0(2N — 2) by the isospin chemical potential 
Hi. The best we can do is to write <j) = (</>o, 0, po, , 0, .., 0) to incorporate both a chiral condensate 
and a charged pion condensate. This will also give rise to off-diagonal entries in the matrix. We 
now set TV = 2 in order to simplify the derivation. Returning to a general N is straightforward. 
We may assume, because of the 0(2) x 0(2) symmetry, that Uk(4>) = Uk(pi,P2) where 



Pi 



P2 = o 



■5. 



If wc denote the fields by double indices 

d 2 U k 



i 

2 ( 

1 2 

Hi 012; ^21; 022 )> we can derive the identity 

d 2 u k 



= ^^SabSij + 4>az4>b] 



idfoj dp b 



dpadpb 



(22) 



(23) 



(24) 



The Fourier transform of R k (x-y) is (2ir) d 8{q+q')R k (q) = Rk,q,q'- If we add bijR k , q , q ' to T ( k J^ - q ql , 
we obtain the matrix 



[rl 2) 



R k ] qq , = (2ir) d 6(q + qi)x 



f d 1 U k + 2p 1 d 2 U k + F k (q) 2^pTp^d 1 d 2 U k 

d 2 U k + F k {q) -2/uqo 

2^pTp2~d l d2U k 2p iqo d 2 U k +2p 2 diU k + F k {q) 

V diUk + Fkfa) J 

where F k (q) = q 2 + R k (q). Using the fact that the inverse in Fourier space satisfies 



F(q 1 ,q') ik F- 1 (q 2 ,q') kj = (2n) d 6( qi - q 2 )S tJ 



,(25) 



(26) 
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we may invert + Rk] , to obtain the full fc-dependent propagator. This may then be inserted 



1 q,a 



into Eq. ([20)1 in order to obtain the RG equation for the effective potential. Going to imaginary 
time we replace the integral over go by a Matsubara sum . This yields 



d k U k = J dkRk 



2(q 2 + R k + d x U k + 2 Pl d 2 U k )(q 2 + R k + 8 2 U k + p 2 d 2 U k ) 

4piP2(did 2 U k ) 2 + (q 2 + R k + d 2 U k )(q 2 +R k + 2 P2 d 2 U k ) + 4/«?w» 
(q 2 + R k + d!U k + 2 Pl d 2 U k ) \(q 2 + R k + d 2 U k + 2 P2 d 2 U k )(q 2 + R k + d 2 U k ) + 4pjuj 2 n 



-4^0,(8^17,0 V + Rk + d 2 U k 



1 



q 2 +R k +d 1 U k 



(27) 



where now q 2 = u 2 +q 2 and ui n = 2itnT denotes the Matsubara frequencies. The difference between 
this and arbitrary N is that there would be 2N — 3 propagators of the form of the last term in the 
above expression. 

Performing the integral over q and returning to arbitrary N, we may simplify this by writing 



d k U k 



ATk d v d 



(d-1) 



d-l \ 

IT 2 - 



27V -3 



NK) 

D{uj 2 ) + ul + fc2 + d x U k 



(28) 



where we have defined 



N(oj 2 ) = 2(k 2 + ujI + dtUk + 2 Pl d 2 U k )(k 2 + uj\ + d 2 U k + p 2 d 2 U k ) - 4p 1 p 2 (a 1 9 2 C/,) 2 

+ (k 2 + + d 2 U k )(k 2 +uo 2 n + d 2 U k + 2 P2 d 2 U k ) + Ap 2 Lo 2 n , (29) 

D(ojI) = (ojI + k 2 + d!U k + 2 Pl d 2 U k ) [(k 2 +uj 2 n + 8 2 U k + 2p 2 d 2 U k )(k 2 +uj 2 n + d 2 U k ) 

+4/i^2] + Ap l p 2 {d l d 2 U k f{k 2 + u 2 n + d 2 U k ) . (30) 

The denominator D(u; 2 ) is a third-degree polynomial in oj 2 . Writing out its roots Pi would take 
several pages. However, they can in principle be found, so we may for the sake of argument assume 
that we have found them. If we denote by Ri the residue of the three poles P{ of l/£>(w 2 ), we may 
write 1 



N(wp A R 



, < + Pi 

1=1 n 



Summing over the Matsubara frequencies yields 



(1 + 2n{y/Pi))Ri 



2\AR 



+ {2N - 3) 



l + 2n(u 3 ) 
2w 3 



(31) 



(32) 



Finding the poles Pi and the residues Ri is done numerically, by use of the residue function in 
Matlab. 
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In the case of zero chemical potential there is only the possibility of a chiral condensate. The 
situation is again 0(2N) symmetric. Eq. (|3"Tj) then reduces to 



N(cul) 1 1 



D{ul) u; 2 l + k 2 + d 1 U k (p 1 ) + 2p 1 d 2 U k (p 1 ) ul + k 2 + diU k (piY 
The equation thus obtained is 



(33) 



where cji = ^Jk 2 + diU k (p\) + 2p 1 d 2 U k {pi) and w 2 = \/k 2 + diU k (pi). 

The large- N limit of Eq. ([32]) is obtained by scaling the fields 4n — > vN <pi . This implies that 
U k NU k on the left-hand side of Eq. (|52"j) . Letting N — ► oo, we obtain 

d k U k = K f- 1 [1 + 2n(w3)] . (35) 



Note that the flow equation (|35|) is independent of the isospin chemical potential. 
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